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ABSTRACT 


A computer program is presented for the analysis of compressible fluid flow across 
a sealing dam, such as is used in a shaft seal. The mathematical model consists of two 
closely spaced parallel rings, one of which is rotating. The analysis is restricted to 
steady, laminar, subsonic, isothermal flow. The effect of rotation on mass flow, pres- 
sure distribution, and other physical parameters is determined. The computer program 
is written in FORTRAN IV. The input and output variables are in English units; printout 
of data in the International System of units is optional. 
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INVESTIGATION OF ISOTHERMAL, COMPRESSIBLE FLOW 
ACROSS A ROTATING SEALING DAM 
II - COMPUTER PROGRAM 

by John Zuk, Patricia J. Smith, and Lawrence P. Ludwig 
Lewis Research Center 

SUMMARY 

A computer program is presented for steady, laminar, subsonic, isothermal, 
compressible -fluid flow across a sealing dam of a radial shaft seal. Output variables 
include mass flow, volume flow, exit Mach number, rotational flow Reynolds number, 
pressure flow Reynolds number, Knudsen number, center of pressure, sealing dam 
force, pressure distribution, power loss, and approximate temperature rise of the gas 
leaking through the seal. The computer program is written in FORTRAN IV for the com- 
puter at Lewis Research Center, which is an IBM 709411/7044 or 7040 direct couple 
computer under IBSYS version 13 using ALTIO. All input, output, and calculations are 
in English units; printout of data in the International System of units is optional. The 
program flow charts are presented in detail and a sample problem is given. 


INTRODUCTION 

In a companion paper (ref. 1) , the analysis is given for the mass flow and pressure 
distribution for a parallel set of surfaces typical of that in the sealing dam of shaft face 
seals (ref. 2). The analysis is for steady, laminar, subsonic, isothermal, compressible- 
fluid flow with rotation of one of the sealing dam surfaces. Mass flow and pressure dis- 
tribution with rotation are compared with the hydrostatic case. Reference 1 shows that 
rotation has significant effects when the speed is hi^ or the radial pressure differential 
across the sealing dam is very small and the speed is moderate. When the sealing dam 
radius is much greater than the sealing dam width, which, in turn, is much greater than 
the film thickness, and the radial pressure flow is sufficiently large, a hydrostatic 
analysis is shown to be valid. Further, in the analysis presented in reference 1, the 
region of validity of the analysis is defined as follows: 



(1) The analysis is valid until the circumferential shear flow causes transition to 
turbulence. 

(2) The analysis is valid until the radial pressure flow approaches a Mach number of 
where y is the ratio of specific heats or until the pressure flow becomes turbulent. 

The objective of this report is to present the computer program for the analysis 
given in reference 1. Input variables include the dimensions of the seal, pressure bound- 
ary conditions, and molecular weight and physical properties of the gas. For specified 
film thicknesses, which are used as parametric input, the output includes mass flow 
rate, pressure distribution, velocity distribution, Mach number, force, center of pres- 
sure, rotational flow Reynolds number, pressure flow Reynolds number, power loss, 
torque, and approximate temperature rise due to the viscous shearing. The computer 
program is written in Fortran IV for the computer at Lewis Research Center, which is 
an IBM 709411/7044 or 7040 direct couple computer under IBSYS version 13 using 
ALTIO. All input, calculations, and output are in English units; printout of data in the 
International System (SI) of units is optional. 


COMPUTER ANALYSIS 
Sealing Dam Model 

A typical face seal is shown in figure 1. It consists of a nonrotating nosepiece held 
either in sliding contact against or in close proximity to a rotating seal seat; the dynamic 
leakage gap (sealing dam) between the nonrotating and rotating surfaces is exaggerated in 



2 



II 


the diagram. Since typical values of the sealing gap are less than 0.0010 inch 
(0.0025 cm), the gap is, in effect, a long narrow slot, and viscous forces are predomi- 
nant in leakage analysis. The model that approximates this sealing dam is shown in 
figure 2. It consists of two parallel, concentric circular rings, separated by a very 
narrow gap, with one moving at a constant speed. A pressure differential exists between 




the inner and outer diameter of the rings, and the fluid is stagnant in the inner and 
outer cavities that boimd the sealing dam. The model formulation is based on the fol- 
lowing restrictions: 

(1) Homogeneous, compressible, viscous, and Newtonian fluid 

(2) Steady and laminar flow with negligible inertia forces 

2 

(3) Bulk modulus, X = (Symbols are defined in appendix A. ) 

(4) Perfect gas fluid 

(5) Negligible entrance region effects 

(6) Isothermal fluid film 

(7) Entrance Mach number, close to zero 


Sealing Dam Analysis 

The governing flow equations for a compressible fluid with constant viscosity in 
vector notation are as follows (ref. 3): 
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Conservation of mass 


+ pv • V = 0 
Dt 


Conservation of momentum (Navier -Stokes equation) 

p 5Y = VP - pVx(VxV) + (A + 2p)V(V • V) + F 
Dt 

Equation of state for an isothermal process 

P = P(P) 

In reference 1 the governing flow equations are applied to the sealing dam model, 
simplified, and solved by 

(1) Utilizing a cylindrical coordinate system (see fig. 2) 

(2) Simplifying the governing equations by using the restrictions listed for the 

sealing dam model 

(3) Applying the appropriate boundary conditions. 

The resulting equations are shown in table I in the form used in the computer program. 

The equations of table I are in the English system of units, and table n shows the 
conversions to SI units. The computer calculations are carried out in the English sys- 
tem of units, but the program contains provision for printout in SI units. 


4 



I 


TABLE I. - FORM OF PERTINENT EQUATIONS FOUND IN COMPUTER PROGRAM 
(IN SAME SEQUENCE AS EQUATIONS APPEAR IN PROGRAM) 



P „inU„ in. /ft) 

Re(P) = = Re, 

(32. 174 <^{T + 460) 

V lbf-sec2/ 


2h 


Q = 13.083 M, scfm 


Power 


jLt AV / 1 h]^ 


\12 HL / 

ft/ 


(42 42MzEH} (Power) 

hp ___ o^ 


MC„ 


Htot = 42.42 (Power), Btu/min 


-'0 


(P - Pmin> “• 


X = Xj + nAx, in. 


P(x) - J 


exp 


[4ta] 




2R^ 


^*^2-^1 
- •'O 




dX 


f 


exp 


X + R^ 


1/2 


2 Elf'll 

dx + pf expl ^ — ' 

\2R^J 


2o2 


where K, = - 


3R|J2‘ 

5(RT 


|Pl - P2KR2 - Ri)L 


(R2 - «!•) 

Torque = (Power) /12 

N \ ft / 
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TABLE n. - CONVERSION TO INTERNATIONAL UNITS 


Quantity 

To convert 

Multiply by 


From 

To 


Length 

in. 

m 

2. 54x10'^ 

Mass 

Ibm 

kg 

0.45359239 

Time 

min 

sec 

60 

Force 

Ibf 

N 

4. 4482216 

Torque 

Ibf-ft 

m-N 

1. 3558179 

Energy 

Btu 

J 

1.05587x10^ 

Power 

hp 

W 

7.4569987x10^ 

Density 

Ibf-secVft'^ 

kg/m^ 

5. 15379X10^ 

Viscosity 

Ibf-sec/ft^ 

N-sec/m^ 

47. 880258 

Gas constant 

ft -Ibf /Ibm -°R 

J/kg-K 

5. 38095 

Specific heat 

Btu/lbm-^R 

J/kg-K 

4. 1865783x10^ 

Temperature 

Of 

K 

K = 5(F + 460)/9 

Temperature difference 

Of 

K 

K := 5F/9 

Pressure 

psi 

N/m^ 

6. 8947572x10^ 

Velocity 

rpm 

rps 

1/60 


ft/sec 

m/sec 

0. 3048 

Area 

in.2 

m2 

m 

6.4516X10"^ 

Mass flow rate 

Ibm /min 

kg/sec 

7. 55987X10"^ 

Volume flow rate at 

scfm 

m Vsec 

4.7194744x10"^ 

standard conditions 




Heat flow rate 

Btu/min 

W 

17. 59783 


Also included are lists of the input variables and program variables along 
with their descriptions and English units, detailed descriptions of the subprograms, 
and a flow chart of the main program (fig. 3). (Listings of the programs are given 
in appendix B; the sample problem, with output options, appears in appendix C.) 


Description of Program 

RSEAL - Main Program . - The main program, RSEAL, performs the primary flow 
analysis. Subroutines are used for secondary operations such as numerical integrations 
and plotting data. Input to RSEAL is by punched cards in the following order: 

(1) Title card: alphanumeric identification of the data (format 12A6) 

(2) NJ card: number of film thicknesses to be analyzed in one running of the pro- 

gram (format 13) 

(3) Film thickness, h^^, cards: 6 per card (format 6F12. 0) 

(4) Data cards: seal dimensions, operating conditions, physical properties of gas 

and logical variables (read by NAMELIST/RINPT/) 

Data are read by NAMELIST to minimize the number of cards required to run a 
second case with the same title and h^^ cards. Input variables are set of zero initially. 
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Consequently, variables that are not changed during the reading of data will be calculated 
by the program. (See the list of input variables for the significance of each variable and 
any restrictions on them.) 

Output data are printed in English units in the following order of groups: 

(1) Program identification: compressible sealing dam with rotation and parallel 

surfaces 

(2) Data identification: printout of title card 

(3) Input data as it appears on RINPT cards 

(4) Calculated constants: gas constant, reference density, speed of sound, length 

of seal, area of seal 

(5) Parameters that vary with h^: mass flow rate, volume flow rate at standard 

conditions, maximum Mach number, pressure flow Reynolds number, rotational 
flow Reynolds number, sealing dam force, center of pressure, dimensionless 
sealing dam force, dimensionless center of pressure, and Knudsen number, 
(Printout of dimensionless sealing dam force and center pressure may be 
skipped by setting RSKIP to TRUE. ) 

(6) Parameters associated with power dissipation: power, total shear heat, appar- 

ent temperature rise, and torque. (Printout of all data associated with power 
dissipation may be skipped by setting TSKIP to TRUE, ) 

(7) Distribution of parameters across the sealing dam (one group for each h^^) : 

distance across seal, pressure, pressure ratio, average velocity, and Mach 
number. 

Printout of data in English units is followed by plots of various parameters. Plots 
appear in standard form with minimum X and minimum Y in the lower left-hand 
corner. Legends at the bottom of the plots give conversion factors for SI units. The 
plots are as follows: 

(1) Power versus h^^ for N ^ 0 

(2) Pressure ratio versus distance across seal for first case for which the model is 
valid. Plot 2 may be suppressed by setting ASKIP to TRUE. 

Following the plots, the data are printed in SI units in the same order as for English 
units. This printout may be skipped entirely by setting NOUI to TRUE. 

RSEAL is divided roughly into eight sections. The first section reads data and 
calculates program constants (cards 41 to 56), 

Section two (cards 82 to 87) calculates the constants needed in the pressure equation. 

Section three (cards 93 to 105) tests input variables. If they are zero, new values 
for them are calculated. Since SPEED and CAPV both represent rotational velocity, 
they must be consistent. If SPEED is read as nonzero, CAPV is calculated from SPEED. 
If SPEED is read as zero, CAPV must be examined. In the case that CAPV is nonzero, 
SPEED is calculated from CAPV. If both SPEED and CAPV are zero, the system con- 
sidered is static. 
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Section four (cards 111 to 139) is the first part of a loop which is done for each h^^j. 
This section calculates a starting value for X, mass flow rate, volume flow rate at 
standard pressure condition, pressure flow Reynolds number, rotational flow Reynolds 
number, maximum Mach number, and Knudsen number. If the maximum Mach number 
is greater than l/V^, this analysis is no longer valid, IHTAG(J) is set equal to 1 as a 
trigger, and no further calculations are made. If the model remains valid, the program 
calculates power, temperature rise due to power dissipation, total shear heat, sealing 
dam force, and center of pressure. The integrations in the force and center of pressure 
equations are done numerically by Simpson's Rule. 

Section five (cards 144 to 154) is the rest of the loop started in section three. 

Section five calculates pressure, pressure ratio, average velocity, and Mach number 
for several points across the sealing dam . 

Section six (cards 158 to 200) writes data in English units. Section seven (cards 204 
to 233) plots the various parameters. And section eight (cards 237 to 295) writes data 
in SI units. 

Numerical constants in the program are for units conversion. 

Subprograms . - RSEAL uses eight subprograms whose listings appear in appendix B. 
These are SIMPSl, SIMPS2, PX, PXX, PRESS, ROT, EXPK, and ARRNG. RSEAL also 
uses two other subroutines that are not standard in IBSYS: SORTXY and PLOTXY. 

SORTXY sorts two numerical arrays. A statement such as CALL SORTXY (X, Y, N) 
results in the array X being rearranged such that X(l) 5 X(2) s . . . 5 X(N) with the 

Y array entries reordered to preserve the (X, Y) pairs. N is the number of elements 
in the X and Y arrays. 

PLOTXY plots two numerical arrays. A statement such as CALL PLOTXY (X, Y, 
KODE, P) produces a printer plot of X versus Y with the minimum X and the mini- 
mum Y in the upper left corner of the page with X increasing down the page and 

Y increasing across the page. KODE tells vtiich plotting options are to be used. For 
example, KODE = 6 gives a plot with most of the grid lines suppressed, with asterisks 
as the plotting character, and with the X and Y scales computed by the plotting routine. 
The array P contains information needed by the plotting routine such as the number of 
points to be plotted, the X and Y scales if the programmer computes them, and 

the frequency of grid lines in the X and Y directions. 

Special format statements are used to print plot titles and plot legends. A pair of 
statements such as 


WRITE (6, 1) 

1 FORMAT (2HPT, lOHPLOT TITLE) 


will print the title 
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PLOT TITLE 


above the plot. A pair of statements such as 

WRITE (6,2) 

2 FORMAT (2HPL, IIHPLOT LEGEND) 


will print the legend 


PLOT LEGEND 


immediately below the plot. 

Listings of SORTXY AND PLOTXY along with details for their use can be found in 
reference 4 or can be obtained from the Instrument and Computing Division of the Lewis 
Research Center. 

SIMPSl is a function subprogram to perform a numerical integration by Simpson’s 
Rule. The integrand is evaluated at interior points by the external function named in the 
calling vector. A statement such as F = SIMPSl (XO, XF, G, K) gives F as the 
definite integral 



G(X) dX 


The integrand is evaluated at interior points by the function G named in the calling 
vector. The interval of integration XO to XF is not divided uniformly. More sub- 
divisions are made where values of the integrand are changing rapidly. If two succes- 
sive evaluations of the integral on a particular subinterval differ by more than 3x10 ^x 
value of the integrand, the subinterval is divided into two subintervals and the integration 
is repeated. If the integration required more than 200 subintervals, the integer K is 
raised by 1 to indicate that the value of the integral is incorrect. 

SIMPS2 performs essentially the same operation as SIMPSl, but it permits carry- 
ing a constant parameter into the integration by means of the calling vector. A state- 
ment such as F = SIMPS2 (A, XO, XF, G, KK) gives F as the indefinite integral 

/•XF 

/ G(A, X) dX 

•'xo 

Use of SIMPS2 permits evaluating double integrals because the external function G can 
call SIMPSl to evaluate the inner integral. KK performs the same function as K, but 
KK is raised by 2 to indicate inaccuracies in the integration. Subprograms SIMPSl and 
SIMPS2 are in general use at the Lewis Research Center. 


9 



PRESS is a function subprogram to evaluate the pressure at any distance from the 
inner radius of the seal. This distance appears in the calling vector. The pressure 
differential equation is solved analytically and the resulting formula is used in the 
program. 



where PX is an external function to evaluate the integrand in the integral 


F 

L 


/ 


Ro *Ri 


(P - 




Similarly, PXX is an external function to evaluate the integrand in the integral 


rRo-Ri 


(P - Pmin)'' ^ 


where P . is the smaller of P^ and P9. 

min 1 ^ 

EXPK is a function subprogram to evaluate exp(K2xy^|^2R2^ where x appears in 
the calling vector. 

ROT is an external function to evaluate the integrand in the integral 
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ARRNG is a subroutine subprogram which arranges two arrays X and Y for plot- 
ting. The subroutine first eliminates cases for which the model does not remain valid. 
It then sorts the remaining data in order of ascending X. It then inverts both arrays 
and interchanges them. The data are now in the form of ordered pairs (Y(N), X(N)), 
(Y(N - 1), X(N - 1), . . ., (Y(l), X(l)) where Y(N) > Y(N - 1) ^ ^ Y(l). Data 

arranged in this form permit plotting with Y as the independent variable in descending 
order. Consequently, the plots appear with the minimum X and minimum Y in the 
lower left comer. 


Input Variables 


The following table lists variables input to the program and their units. Arrays are 
given with their dimensions. 


Variable 

Unit 

Description 

TITLE (12) 


Data identification 

NJ 


Number of mean film thicknesses (NJ 50) 

HMEAN(J),-J=1,-NJ 

in. 

Mean film thicknesses 

L 

in. 

Width or mean circumference of sealing dam 
of flow 

SPEED 

rpm 

Rotational speed. Either SPEED or CAPV 
may be set. If both are set, calculate CAPV 
from SPEED 

CAPV 

ft/sec 

Surface speed. 

MOLWT 

lbm/(lb~mole) 

Molecular weight of the gas 

PI 

psi 

Pressure at inner radius of seal 

P2 

psi 

Pressure at outer radius of seal 

T 

Op 

Isothermal reference temperature 
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Variable 


Unit 


Description 


R1 

R2 

RHORO 


RHORF 


MU 

CP 

GAMMA 

NGRID 

ASKIP 


RSKIP 


TSKIP 


NOUI 


in. 

in. 

Ibf-sec^ft^ 

Ibf-secVft'^ 


Ibf -sec/ft ^ 
Btu/lbm-°R 


Inner radius of seal 
Outer radius of seal 

Reference density at inner radius of the seal. 
If RHORO is read as zero, program cal- 
culates RHORO. 

Density at mean radius used in calculating 
the rotational Reynolds number. If RHORF 
is read as zero, the program calculates 
RHORF. 

Absolute viscosity of gas 
Specific heat of gas 
Ratio of specific heats 
Number of steps across seal (max = 20) 
Logical variable. If ASKIP = TRUE, pro- 
gram skips calculation and printout of x, 
pressure, average velocity, Mach number, 
and pressure ratio. It also skips plotting 
of X versus pressure ratio. 

Logical variable. If RSKIP = TRUE, program 
skips calculation and printout of dimension- 
less center of pressure and force. It also 
skips plotting of dimensionless F. 

Logical variable. If TSKIP = TRUE, program 
skips calculation, printout, and plotting of 
power. 

Logical variable. If NOUI = TRUE, program 
makes no conversion to SI xmits. 


Program Variables 

The following table lists the variables used in the program in the approximate order 
of their appearance, any limitations on the variables, and their units. Variables marked 
with an asterisk are printed as output data. 

PI n = 3 . 1415927 

RUNTV ft-lbf/(lb-mole)°R Universal gas constant 
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Variable 


Unit 


Description 


ZERO 


PREF 

psi 

RRl 

in. 

MCUT 


NN 


JMOD 


RDIF 

in. 

PDIF 

psi 

AREA* 

in.^ 

R* 

lbf-ft/lbm-°R 

RHOl* 

Ibf-sec^/ft^ 

A* 

ft/ sec 

RHOREF* 

Ibf-sec^/ft^ 

OMEGA 

rad/ min 

ECONST 

1/in.^ 

EA 

Ibf^/in.^ 


AA Ibf^/in.^ 


Input variables are set to 0. 

Smaller of PI and P2 in the COMMON 
block 

Value of R1 in the COMMON block 
1/ = upper limit for which model is valid 

Number of grid points (1 ^ NN < 21) 

Integer number of pressure distributions that 
will fit on one page. 

Distance across seal = R2-R1 
Total pressure drop across seal 
Face surface area of seal 
Gas constant 

Density at inner radius of seal 
Speed of soimd 

Calculated density at mean radius of seal 
Rotational velocity 

o 

-Kj/ 2R| (where Kj is defined in appendix A) 


2 ^1 
p| exp 




DELX 

x(j, D* 

MDOT(J)* 

VAV 

MACHMX(J)* 

REP(J)* 


in. 

in. 

Ibm/ min 
ft/sec 


Distance between two successive grid points 
Distance from inner radius of seal 
Mass flow rate 

Average velocity at outer radius of seal 
Maximum Mach number for given HMEAN 
Pressure flow Reynolds number 
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Variable 


Unit 


Description 


RER(J)* 


Rotational flow Reynolds number 

Q(J)* 

scfm 

Volume flow rate at standard conditions 

KN(50)* 


Knudsen number 

IHTAG(J) 


Numerical flag 

IHTAG(J) = 0 implies model remains valid 
IHTAG(J) ^ 0 implies model is invalid 

POWER(J)* 

hp 

Power dissipated by viscous shearing 

DELTJ(J)* 

°F 

Apparent temperature rise due to power 
dissipation 

TORQUE (J)* 

Ibf-ft 

Torque 

HTOTAL(J)* 

Btu/min 

Total shear heat of system 

F(J)* 

Ibf 

Equivalent force 

XC(J)* 

in. 

Center of pressure 

K, KK 


Numerical flags that indicate whether the 
numerical integrations in the calculation of 
F(J) and XC(J) are accurate 

P(J, I)* 

psi 

Pressure at X(J, I) 

PRAT(J, I)* 


Ratio of P to P^. at X(J, I) 
min 

UAVRG(J, I)* 

ft/sec 

Average velocity at X(J, I) 

MACH(J, I)* 


Mach number at X(J, I) 

FBAR(J)* 


Dimensionless force 

XCBAR(J)* 


Dimensionless center of pressure 

JJ 


Counter for valid cases. If JJ = 1 
modulo JMOD, the printer will skip to 
a new page. 

XPT(50), YPT(50), 


Utility arrays used in sorting and plotting 

XPLOT(50), YPLOT(50) 

1 

data 

PP(61) 


Array needed by plotting subroutine - see 
description of PLOTXY for details 

KOBE 


Plotting code - see description of PLOTXY 
for details 

NP 


Number of points in a plot (1 < NP ^ 50) 

I 


X index 

J 


HMEAN index 
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Subroutines 


Name 


Call Common Program 

vector block variables 

variables variables 


Description 


PRESS - pressure function 


ROT - integrand in integral 



exp 


Kj^X + R^ 


2 V Ro 


-dx 


X + R. 


PX " integrand in integral 
'R2-R1 
'0 


r 


(P - P . ) dx 


R1 

X 


PREF 

RRl 

EA 

AA 


PREF 
A, B, 
C, D 


R1 


QQ 


K 


Q 


Distance from inner radius of seal 
Reference pressure 
Inner radius of seal 


p; exp —ti 

j- O V T5 


2 \R. 


2 Ki 2 

exp - Pf exp -i — ^ 

2 2 ^ 2 \R2/ 

^ 2 



dx 


X + R. 


Inner radius of seal 


K. /X + R. 
exp — ' 



dx 


X+ R, 


Numerical flag indicates whether 
numerical integration is 
accurate 
p2(x) 

Inner radius of seal 

Distance from inner radius of seal 

Radial distance = R1 + x 


Distance from inner radius of seal 

Reference pressure 

Dummy variables to fill COMMON 

block 

P - P . 
min 
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Name 

Call 

Common 

Program 


vector 

block 

variables 


variables 

variables 



Description 


PXX - integrand in integral 



-Ri 




X 


EXPK - special exponential X 

function 


ARRNG - arranges arrays in X 

proper order for plotting Y 

XP 

YP 

N 

I 


PREF 


Distance from inner radius of seal 
Reference pressure 


A, B, 

C, D 

Y 


ECONST 

YY 


Dummy variables to fill COMMON 
block 


(P - Pmin)X 

Square of distance from center of 
seal 


(Kj/2)/r| 


Distance from center of seal 


Y 


T 


exp 



Input array of independent variables 
Input array of dependent variable 
Sorted and inverted array of new 
independent variable 
Sorted and inverted array of new 
dependent variable 
Number of elements in input arrays 
Number of elements in sorted 
arrays 

Temporary storage for sorting 
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Do statements 200 to 400 for each J 
J = 1. NJ 



(a) Initial steps. 




Figure 3. - Flow chart of Main Program. 


lb) Main calculation. 
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(d) Plot routine. 
Figure 3. - Concluded. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, Jime 24, 1969, 
120-27-04-90-22. 
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APPENDIX A 


SYMBOLS 


X. , .22 

A cross-sectional area, in. ; m 

a speed of sound, ft/ sec; m/ sec 


Q net volume flow rate, scfm; 
std. mVsec 


C specific heat at constant pressure, 
^ Btu/(lbm)(°R); J/(kg)(K) 

C specific heat at constant volume, 
Btu/(lbm)(°R); J/(kg)(K) 

pj 

D/Dt material derivative, '^ + ^ 

F sealing dam force, Ibf; N 

F dimensionless force, 

F/(P2 - P^)(R2 - Ri)L 


R mean radius, (Rj + R2)/2, in.; m 

AR sealing dam length, R 2 - Rj, 

in. ; m 

Rej^ Reynolds number in radial direc- 

tion, puh/p 

Re^, Reynolds number in circumferen- 

tial direction, PRS2h/p 

tfr gas constant, ^/m, 

ft-lbf/(lbm)(°R); J/(kg)(K) 

(R universal gas constant, 1545.4 


F body force vector 

h film thickness, nominal, in. ; m 

-(3R|fi2/5(RT) 

L sealing dam width, in. ; m 

M Mach number 

M mass flow, Ibm/min; kg/sec 

AM change in mass flow 

m molecular weight of gas, 

lbm/(lb-mole); kg/kg-mole 

N speed of rotation, rpm 

n integer 

n 

P static pressure, psi; N/m 

AP pressure differential, psi; 

N/m^ 

Pmin smaller pressure of two pressure 
boundary conditions, psi; 

N/m^ 


ft -Ibf / (lb -m ole) (°R) ; 

8. 3143 j/(kg-mole)(K) 

r radial direction coordinate 

T temperature, °F; K 

T average temperature, °F, K 

U pres sure -flow reference velocity, 

ft/sec; m/sec 

u velocity in r -direction or 

x-direction, ft/sec; m/sec 

V reference shear flow velocity, 

ft/sec; m/sec 

V velocity in 0-direction, ft/sec; 

m/sec 

w velocity in y -direction, ft/sec; 

m/sec 

W reference velocity across film 

thickness, U(h/AR), ft/sec; 
m/sec 
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Xp center of pressure in radial or 

x-direction, in. ; m 

dimensionless center of pressure, 

X^/(R2 - R^) 

X coordinate in pressure gradient 

direction 

y coordinate across film thickness 

z shear flow coordinate in Cartesian 

system 

y specific-heat ratio, Cp/C^ 

0 circumferential coordinate 

X second viscosity coefficient or 

coefficient of bulk viscosity 

M absolute or dynamic viscosity, 

(lbf)(sec)/ft^; (N)(sec)/m^ 

2 

V kinematic viscosity, ft /sec; 

m^/sec 


p density, (lbf)(sec^)/ft^; kg/m^ 

Q, angular rotational velocity, 

rad/sec 

V Del operator, i + -J- ^ j + ^ k 

Subscripts: 
av average 

h based on film thickness 

m mean 

r based on radius 

ref reference 

0 reference 

1 inner radius or inlet 

2 outer radius or exit 
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APPENDIX B 

PROGRAM LISTING 

IbFTC RSEAL DEBUG 

COMPKESSIBLE FLOW SEALING DAM ANALYSIS WITH ROTATION) 


LOGICAL ASK1P,RSK1P,TSKIF,N0UI 
REAL MDOT,MOLwT,MUfMACri,MACHMX,MCUTfKN,L 

DIMENSION XPLOK5 0) ,YPLCT(5 0) ,XPT<50) ,YPT (50) , PPC61 i , ZERO! 171 . 

I UlI30),TlTLE(i2) 

DIMENSION F(50) ,XC(50) , MOOT (50) ,XCBAR(50) ,FBAR( 50 ) , Q( 50 ), 

1 POWER ( 50) .HI (50) .H TOTAL (50) ,D£LTJ(50 ) .TORQUE ( 50) ,HMEAN(50), 10 

2 MACHMX(50i ,KEP(50) tIHTAG(50) ,KN(50) ,RER(50) 11 

DIMENSION X( 50,21) ,P( 50,211 .UAVRG (50 ,21 ) , PRAT ( 50 , 21 ) • MACH( 50, 2 1 ) 12 

COMMON/ INTGRL /PREF ,RR 1 , EA ,AA, S 13 

COMMON/ ECON/ECONST 14 

EXTERNAL PXtPXX.ROT 15 

IMAMELIST/RINPT/L , SPEED, CAPV ,MOLWT ,P1 ,P2,T ,R1, R2, RHORO, 16 

1 RH0KF,MU,CP,GAMMA,NGR10, ASKIP,KSKIP,TSKIP,NOUI 17 

data pi ,KUNIV/3. 1415927,1545.4/ 18 

equivalence (ZER0(1),M0LWT) , (ZER0(6), Pi), ( Z ERO( 1 1) ,L ), 19 

1 (ZERU(2) .SPEED) , (ZER0(7), P2) , ( ZERO( 1 2) , T ), 20 

2 IZ£KU(3) .RHORO) , (ZER0(8), Rl), ( ZERO( 13) ,MU) , 21 

3 (ZEKU(4i ,RHOHF) , (ZER0(9), R2), ( ZERO( 14) , CAP V 1 , 22 

4 (ZEKU(5i ^GAMMA) , ( ZERO ( 10 ) , CP ) 23 

DO 90 1 =i ,14 24 

90 ZEKO( 1 ) = 0. 25 

C 26 

C READ input DA TA , l ALCULA TE PROGRAM CONSTANTS, AND WRITE INPUT 27 

C CONDITIONS 28 

C 29 

C DATA CARDS 30 

C TITLE “ DATA IDENTI F I CA TI ON - 1 CARD (FORMAT 12A6) 31 

C 32 

C NJ - NUMBER OF FILM THICKNESSES (FORMAT 13) 33 

C 34 

C HMEAN - MEAN FILM THICKNESSES - 6 PER CARD (FORMAT 6F12.0) 35 

C 36 

C $KINPT - SEAL DIMENSIONS, OPERATING CONDITIONS, PHYSICAL 37 

C PRUPERTiES OF GAS, LOGICAL VARIABLES 38 

C (READ BY NAMELI ST/RINPT/) 39 

C 40 

READ (5,3) TITLE 41 

READ (5,1) NJ 42 

READ (5,2) (HMEAN( Ji , J = 1 ,NJ) 43 

loo WKl IE ( 6, 10) 44 

read (5.K1NPT) 45 

PREF = AMIN1(P1,P2) 46 

RRl = Kl 47 

MCUr = i./SQRT(GAMMA) 48 

NN = NGKID^'I 49 

JMOD = 59/(4+NN) 50 

AREA = P1*(R2««2-R1**2) 51 

RD1F=R2-K1 52 

PDlF=A[)S(Pi-P2) 53 

R = RUNI V/MOlwT 54 

RhOi= P l/k/(H-460.) <^4,475o636 55 

A= oQRT(GAMMA*R*n>460.)<'32.1 74) 56 

WRITE (6,54) TITLE 57 

WRITE (6,11) P2 ,P1 , I ,MU,MOlw T,GAMM A, R2 ,R1,C, RHORO, RHORF.NN, 58 

1 SPEED, CAPV, CP, ASKi P, RSKIP.TSKIP 59 

IF (NOUI) GO TO 110 60 
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61 

CONVekT INPUT UATA TO I NTEkkATI UNAL UNITS 62 

63 

JKU = ALPHA 64 

UH2J = P 6. 694 757^ t3 65 

UU3J = P l»t3. 894 757263 66 

UI(4i = <T+460.)/1.8 67 

UK 5)= MU447. 860258 68 

UT(6J= MULWT 69 

UK 71 = GAMMA 70 

UK8J = k2^2. 546-2 71 

UK9) = R 1*2. 546-2 72 

UKiO) = L*2.54E-2 73 

UKIU = RhUKa*517.202o 74 

UK 12) = RHURP*5l 7.2026 75 

UK 13) = SPcEO/oO. 76 

UK14)= lAPV*.3048 77 

UK 15) = CP*.4l 865783E4 78 

79 

UEIERMINE RUTATIUN CONSTANTS 80 

81 

110 OMEGA = 2.*PI*SPE6l) 82 

ECUiMST = -.3*OMtGA**2/R/<T+460.)/3600./32.l74/144. 83 

EA = 6XPMR1**2)*PI**2 84 

AA = EXPK(R2**2)*P2**2-£A 85 

K = 0 ■ 86 

S = 8IHPS2(R1 ,0. ,R0IF ,ROT ,K) 87 

IF (K.N6.0) WRITE (6,21) 88 

UE6UG OMEGA, 6CUNST,EA,AA,S 89 

90 

TEST INPUT PARAMETERS ANU DETERMINE OPERATING CONDITONS 91 

92 

120 IF (L.EQ.O.) L=P1*(R1 + R2) 93 

IFIRHURO.NE.O.) RH01=RH0R0 94 

IF (RHUKF.NE. 0. ) GO TO 123 95 

RFUREF = Pk£SS{RUlF/2.)/R/(K-460.)*4. 4756636 96 

GO fO 124 97 

1/C3 RFOREF = KHOkF 98 

124 IF (SPEED. EU.O.) GO TO 121 99 

CAPV = P1*SPEEU*( RH-R2) /720. 100 

GO TO 122 101 

121 IF (CAPV.NE.O.) GO TO 122 102 

TSKIP = .TRUE. 103 

CP = 0. 104 

122 SPEED = CAPV*720./PI/(R1+R2) 105 

106 

DEItRMlNE STARTING VALUES FOR X AND H 107 

DETERMINE X INCRIMENT, M(UOT), Q, POWER, MAX MACH NUMBER, 108 

PRESSURE REYNOLD'S NUMBER, DELTA T, HEAT(IN), AND HEAT(TOTAL) 109 

110 

200 DO 4C0 J=i,NJ 111 

IhTAG(J) = 0 112 

DElX = RulF/FLUAT(NGRID) 113 

X U, 1) = 0. 114 

MDDT(J) = -13.405833*PI*HM£AN(J)**3*RH01*AA/MU/PI/S 115 

DEBUG HMEANI J) ,MOOT( J) 116 

VAV = ABSIMDOTIJ) )*Pi/L/RH01/HMEAN(J)/PREF/ 13.405833 117 

MaCHMX(J) = VAV /A 118 

IF (MACHMX(J) .LT.MCUT) GO TO 223 119 

IhTAGIJ) = 1 120 
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GO TO 400 121 

223 REP(J) = 2. 4PRtf*VAV*HM£AN<J)/MU/R/ (H-460.)/2. 6811667 122 

222 .<£R(J) = RHOkfcF*GAPV*HMfcAN<JI /MU/12. 123 

U(J) = 13.083*MJOT( J) 124 

KN(J) = 2.96«=MACHf-1X( J) /REPi J) 125 

230 IF (TSKIP) GJ fO 240 126 

POwER(J>= MU*AKEa*CAP V** 2/HMeAN( JI/6600. 127 

JELTj(J) = 42.42*P0h(ER( JJ /A8S (MOOT Ul )/ CP 128 

HTOTAL(J) = 42.42*P0wER( J) 129 

rUi<UUElJ)= POWER! J)*3.3E4/i»PEED 130 

C 131 

C DETERMINE FORCE AND CENTER OF PRESSURE 132 

C 133 

240 K = 0 134 

KK = 0 135 

F(JI= SIMPSUO. ,RUIF,PX,KI*L 136 

XC(JI = SIMPSKO. ,ROIF ,PXXfKK)/F( J)*L 137 

IF (K.NE.O) WRI TE (o,24) HMEAN(J) 138 

IF (KK.NE.OJ WRITE (6,25) HMEAN(J) 139 

C 140 

C determine pressure, film thickness, pressure RATIO (P/Pl), 141 

C AVERAde VELOCITy, AND MACH NUMBER AT EACH GRID POINT 142 

C 143 

300 IF (ASKIP) GO TO 400 144 

P (J, 1) = PI 145 

DO 320 1=1, NN 146 

IF ( I .E(J. 1) GO TO 310 147 

XU, I) = FlDATII - li <‘OtLX+ X( J, U 148 

P(J,I) = PR£SS(X( J,I) ) 149 

31DPRAT(J,1) =P(J,i)/PREF 150 

UA\/KG(J,I )= AttS(MDQT( J) )*P1/L/RH01/P( J, I) /HME AN( J )/ 13.405 833 151 

MACH(J,I) = UAVRG(J,I)/A 152 

320 CLiNT INUE 153 

400 CUImTINUE 154 

C 155 

C WRITE RUUTINE 156 

C 157 

600 WRITE (6,1^) RHDREF ,RHU1 , A, SPEED, R,L, AREA, CAPV 158 

WRITE(6,13) 159 

IF (.NOT.RSKIP) WRITt(6,i4) 160 

WRITE (6,55) 161 

IF (KSKIP) GU to 603 162 

WRITE (6,56) 163 

C 164 

603 DU 610 J=l, NJ 165 

WRITE (6,23) HMEAN(J) 166 

IF (IHTAG(J).Ey.O) Go TO 609 167 

WRITE (6,20) 168 

GO TO 610 169 

609 WRITE (6,15) MDOT( J) ,0( J) ,MACHMX( J) ,REP( J ) ,RER( J) ,KN( J) ,F( J ), XCU ) 170 

IF (RSKIP) GU TO 610 171 

C 172 

C NUKMALI/.E F AND XC 173 

C 174 

FBAR ( J)=F ( J)/PDIF /kDIF/L 175 

XCBARIJ )=XC( J)/kDlF 176 

WRI(t(6,16) XCbAR(O) ,FBAK( J) 177 

610 CONTINUE 178 

IF (TSKIP) GU TU 612 179 

WRITE (6,17) 180 
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00 oli J=1,WJ 181 

riKirt (6,23) Hi-ltAlNUJ) 182 

lF{lHrAG(J).E;U.0)G0IU613 18 3 

WKITE(6,20) 184 

GU TO 611 185 

6iJ rtkiTtl6,22) FUkfcfK J) ,HTuTAL(J) ,OtLTJ( J) ,TCR0UE( J) 186 

611 CONTINUE 187 

612 IF (A6K.IP) GO To 1 00 188 

189 

WRITE X AND FKE66UKE 01 STtU BU f IONS 1^0 

191 

jj = 0 192 

DO 620 J=l, NJ 193 

IF (IHIAGIJ).NE.O) GU TO 62 0 194 

JJ = JJ+1 195 

IF (MOO (JJ.JMOO) .Ew. U WRITE (6,57) 196 

WftITE (6,18) HMEAN(J) 197 

WRITE (o,I9) (XU ,I ) ,PKAT ( J,I ) ,P( J, I J ,UAVRG(J , I J , MACH( J , U, 198 

i 1=1, NN) 199 

620 CONTINUE 200 

201 

PLOI ROUTINE 202 

20 3 

TOO K(jut = 6 204 

PP(3) = 0. 205 

PP(4) = 0. 206 

DO 701 J=1,NJ 207 

IF ( IHTAGl J) .ME. U) rl ME AN U ) =-HMEAN 1 J J 208 

701 CONTINUE 209 

720 IF (TSK.IP) GO TO 730 210 

If ( SPEED .Ew. U. ) Go TO 730 211 

CALC ARRNG{HMEAN,PUwER,XPLOI,yPLQT,NJ,NPJ 212 

PP( 1) = NP 213 

WRITE (6,34) 214 

CALL PLOTXY(XPLOT,yPLOT,KODE,PPJ 215 

write (6,43) 216 

730 IF (ASRIP) Go TO 750 217 

DU 741 J=1,NJ 218 

IF ( IHTAu ( J). NE. 0) GO TO 741 219 

WRITE (6,37) HMEAN(J) 220 

DO 740 1=1, NN 221 

XPK 1 ) = X( J,i) 222 

yPT( I) = PRAT(J,I ) 223 

740 oONl INUE 224 

CALL ARKNG ( XP I , 7 PT , XPL OT , Y FL UT ,NN, NP ) 225 

PP(l) = NP 226 

lALl PLUTXY(XPlOT,YPLOT,K,ODE,PP) 227 

WRITE ( 6 , 46 ) 228 

GO TO 750 229 

741 CONTINUE 230 

750 DO 751 0=1, NJ 231 

IF ( IHTAGl J).i')E.O) HMEANl J ) = -HMEAN( J) 232 

751 CONTINUE 233 

234 

CONVERT 10 INTERNATIONAL UNITS AND PRINT 235 

236 

800 IF (NOUI) GO TO lOU 237 

UI(21)= RH0KtF*517.2026 238 

UI(22)= RHU1*51 7. 2026 239 

U I( 23) = A*. 3 04 8 240 
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UKi*.) = SPEEl)/60. 241 

uKzaJ = k#5.3au95 242 

UH26)= L«2.54E-2 243 

UU27I= AKEA<‘6.431oh-4 244 

UH2a)= uAPV'«=.304a 245 

WRITE I6,1(J) 246 

WRITE I 6f 54i title 247 

WRITE (6.30J (UI m tl=2>12) .NN.IUim ,I=13.I5),ASKlP,RSKIPt 248 

1 TSKIP, (Uiai ,1=21,281 249 

WRITE 16,13) 250 

IF I .NJ T.RSKIP) WRI TE (6,14)‘ 251 

WRITE (6,51) 252 

IF I.NOT.RSKIPJ WRITE 16,56) 253 

00 aiO J=1,NJ 254 

UF = HMEAM J)«=2. 54E-2 255 

WRITE 16,23) uH 256 

IF < IHTAGlJ).EQ.O) GO TO 809 257 

WRIIE (6,20) 258 

GO Tu aiO 259 

a09 UM = MOOT ( J)<=. 755 98 7E-2 260 

00 = oIJ)*. 4 7194 744E-3 261 

OF = F( J)*4. 4432216 262 

UX(, = XC I J )*2. 5tE-2 263 

WRITE (6,13) UM,Uw,MACHRX(J) ,REP(J) ,RER(J),KN( J),UF,UXC 264 

IF( .NOT.kSKIP) WRITE(6,16) XCbAR I J) ,F 8AR( J ) 265 

alO lOiNTINOE 266 

IF (TSKIP) GO TO al2 267 

WHITE (6,52) 268 

00 dll J=1,0J 269 

OF = HMEAfiil 34E-2 270 

WRITE (0,23) UH 271 

lF(IriTAG( J).fM£.U) 30 TO 813 272 

OP = Pu«ER( J)*7.45699d7E2 273 

OFT = rlIUTAL(J)*l 7. 397833 274 

G( = OcLTj(J)/1. 8 275 

0 IKK = rOHQUE(J)>i‘l. 3558179 276' 

WKITE(o,^2) UP,GHT,UT,OTRK 277 

GO TO dll 278 

613 WRITe(o,20) 279 

811 CuhT Ii>tUE 280 

812 IF (ASKIP) GO TO 100 281 

JJ = 0 282 

OC 821 J=1,HJ 283 

IF ( IHTAG( J) . i\e. 0) GO TO 821 284 

Jd = JJ+1 285 

IF (rtuO ( JJ,JMLO) .EU. 1) WRITE (6 ,57) 286 

OF = HnEAN( J)#2. 34E-2 287 

wRIIc (6,33) UH 288 

00 820 1=1, NN 289 

0/ = X( J,I)*2. 54E-2 290 

UP = P( J, I )*6. 894 73 72E3 291 

0V= UAVKG ( 0,1 )<=. i048 292 

WRITE (6,19) UX,PHAT( 0,1) ,UP,UV,MACH( 0,1) 293 

820 CONTINUE 294 

821 CONTINUE 295 

GO TO 100 296 

C 297 

1 FORMAT (13) 298 

2 FORMAT (6F12.6) 299 

3 FORMAT (12A6) 300 
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10 FORMAT ( Irii, 6U^oOrtPKfcSSl 6Lfc DEALING 0AM WITH ROTATION ANO PARALLEL 301 

1 SUKFA^ESJ 302 

11 FORMAT ( iHUflZHl NPUT DATA - ,/ t 303 

1. ihOf oX , 7HP2 ,PS1 A ,IUX, /riPl I Pil A ,iOX,7HT ,UElj F , lOX , 30A 

2 2UH VI jLOGi I Y ,Li3- itC/FT2 f lOX, 16HM0LECUL AR WEI GHT, ISX , 5HCP/CVi /i 305 

J IH , 2G i 7. :?,F 13.0 ,G23. 5 ,G3I .5 fG28.5 ,/ ,1H0 ,4X, 9HK2, INCHESi 8X, 306 

4 9ri« i ,I NGHtS, 9X,8HL ,1 NChE S, 10X,15HKH0,LB-SEG2/FT4 f iOX , ’ 30 7 

5 2uHKHU(K0n ,Ld-3tL2/FT4,10X,l7HNG OF GRID POINTS, /,1H ,3G17.5, 308 

6 G21 .5, Ga7. 5, 1 22 , IH0,6X ,5HN ,KPM,11 X ,8HV ,FT7 SEC, 9X, 309 

7 15HLP, Bl 0/Lb-uEG K,I3X,6HSMP A,6X,6H5K,iP R,8X,6HSKIP T,/,LH , 310 

0 2G1 7.5,G18.3,l19,2L12) 311 

12 format ( IHO, 1 TridEGiN OUTPUT OAT A, / ,1 HO ,7 X , 312 

1 20tlKH0(KUn,LD-5EC2/FT4,4X,18HRH0m ,L3-SEC2/FT4,6X, 313 

2 2irtA(oUUNU ;>PEt J ) ,F r/StC ,6X,9HSPEED ,RPM,/ , IH , 2G22 .5 , G26 . 5, G23 .5, 314 

j /,1H0, 3X,3UHGAb CONi TAivT ,F T- LB/ LB ( M) -UEG R, lOX , 8 HL , I NCHES , 12 X, 315 

4 EHAREA ,iN2,12X,3HV, FT/SEC, / ,IH ,G24 . 5 , G3 1 . 5, 2G20 . 5 ) 316 

13 format ( 1HU,2X, 9HMEAN F IL M ,4X ,6 HM ( UOT I ,6 X , IHU . 8X , 4HMACH , 5X , 317 

1 5HRE IP ) , bX,5HkE I RJ ,5X, 7HKNUDSEN, 7X ,1HF,9X, 2HXC) 318 

14 FORMAT I 1H+, 100X,2HXC ,bX,lHF) 319 

15 format I lH+,iiX,2Gll. 3,2F9. 3,2F10.3,F12.3 ,G11 .3) 320 

16 FUrMaT (1H+,93X,^F10.3) 321 

17 format (lH0,l6HMtAN F IL M, 1 NCH ES ,8 X , 10 HPOW ER , H . P. , 6X , 322 

1 13HSHEAR HEAT, BTU/MIN, 5X, 12HUEL(T» , DEG F ,8X, 12HTQRQU E, FT-L B) 323’ 

18 Format IIHO.IIHMeAN film = ,GI l .3 ,2X,6HINCHES,/ , IHO, 5 X , SHX, inches, 324 

1 4X,8HP/P(M1NJ,7X,5HP ,PSI ,6X,12HU(AV) , FT/SEC, 4X,7HMACH NO) 325 

19 FORMAT (IH ,0 14.3,40 14.5) 326 

20 format (1H+, 77X, 18HANALYSI S NUT VALID) 327 

21 FORMAT I 1H0,39HNUMErICAL INTtoRATIUN IN S IS INCORRECT) 328 

22 FORMAT (1H+, 15X,4o2U.3) 329 

25 format (IH ,G11.3) 330 

24 format (IH0,26HF0R MEAN FILM IHICKENESS =,G15.5,3X, 331 

1 ItHF IS INCORRECT) 332 

25 format (1H0,26HFUR MEAN FILM THICKENESS =,G15.5,3X, 333 

1 15HXC IS INCORRECT) 334 

30 FORMAT (2HPT,55HPL0r OF M(OCT) AND J VS H(MEAN) WHERE U = 13.083*M 335 

llOoT) ) 336 

32 FORMAT I 2HPT,23HPLOT OF X(C) VS H(MEAN)) 337 

33 FURMAT ( 2HP T , 28HP LU T OF MDOT(BAR) VS H(MEAN)) 338 

34 format I 2HP T, 24HPL0T OF PUwER VS HI MEAN)) 339 

37 FORMAT (2HPT,35HPLOT UF P/PIMIN) VS X FOR HIMEAN) =,G15.5) 340 

36 FORMAT (2hPT,25rlPLUT UF FII3AR) Vi HIMEAN) ) 341 

40 FORMAT (2HPL,50HXC ON VERTICAL SCALE - HIMEAN) UN HORIZONTAL SCALE 342 

1 ,/,2HPL,/,2HPL,10X,5oiIXC IN INCHES - 10 CONVERT TO METERS, MULTIP 343 

2LY oY 2.54E-2 , / ,2HP L , / , 2H PL , lOX ,6 1 HHI ME AN ) IN INCHES - TO CQNVER 344 

3T TO HETERi, Multiply by 2.54E-2 ) 345 

41 Format I2HPL, 34HMIbAR) On vertical SCALE - HIMEAN) ON HORIZONTAL S 346 

ICALc , /, 2HPL,/,2MPL, 10X,23HMll3AR) IS 0 IME NS I UNL ES S , / , 2HPL, / , 2HPL , 347 

2 loX,61HriIMEAN) IN INCHEi - TO CONVERT TO METERS, MULTIPLY BY 2.54 348 

3E-2) 349 

42 FORMAT 12HPL,54HF (BAR) ON VERTICAL SCALE - HIMEAN) ON HORIZONTAL S 350 

1CAlE,/, 2HPL,/,2HPL,1UX,23HFIBAR) IS U IMENS I ONL ES S , / , 2HPL, / , 2HPc , 351 

2 iOX.blHHIMEAN) IN INCHES - TO CONVERT TO METERS, MULTIPLY BY 2.54 352 

3E-2) 353 

43 format I 2HPL, 53HP0wER UN VERTICAL SCALE - HIMEAN) ON HORIZONTAL SC 354 

lALE, /, 2HPL,/ , 2HPL ,10X .olHPOWER IN HORSE POWER - TO CONVERT TO WATT 355 

2i, MULTIPLY BY 74 5. 7 ,/ ,2H PL ,/ ,2HPL , 10 X, 6 1 HH I ME AN) IN, INCHES - TO C 356 

BUNVErT to meters, multiply by 2.54E-2,/,2HPL,/,2HPL,50HF0R SHEAR H 357 

4EAT IN bTU/MIN, MULTIPLY POWER BY 42.42) 358 

46 FORMAT I2HPL , 50H P /P I MI N) ON VERTICAL SCALE - X ON HORIZONTAL SCALE 359 

1,/,2HPL,/,2HPL,10X,25HP/P(M1N) IS 0 1 MENS I ONlES S , / , 2HPL , /, 2HP. , 1 OX, 360 
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l fiaHX IN If4CHcS - TO CONVEKT TO METERSi MULTIPLY BY 2.54E-2I 361 

i>U FORMAT (IH0,3X,7MP2,N/M2, 362 

1 SX, 7HPi,N/M2,10X,7rlT,Ut0 K ,1 OX , l&HV I SCOS I T Y, N-S/ M2t I OX , 363 

2 loHMOLbCULA^ WE I OR T , 14X, 5HCP/C Y , / , IH ,G 14. 5, G1 7. 5, F 14.0, G24. 5, 364 

3 G2o.3,G2 7. 3,/,H U,4X ,9HR2 , ME TEKS ,7X,9HR1 .MET ERS , 9X , 8HL.M ETERS, 365 

4 13X,9HKH0,KG/M3,12X,14HRH0(RUT) ,KG/M3,10X, 366 

5 17HN0 OF GRIO POINTO./.IH ,G14 .5 ,G17 .5 ,G Id .5 , G2l . 5, G22 .5 , I 22, /, 367 

6 lMU,dX ,5riN,KP0,12X,5HV,M/S,9X,i3HCP,J/KG-0£G K, 368 

7 16X,6MSKIP A,6X,6riSKIP K,6X,6HSKIP T,/,1H ,G 15. 5 , G1 7.5, 369 

8 Gid.5,i.20,2Ll2,/,lHU,17rib£GIN OUTPUT UAT A ,/ , 1 HO , 5X , 370 

9 i4HKHU(ROTi,KG/M3,iOX,l2HRHOm ,KG/M3,10X, 371 

X lariAIbUUNO SPEEu) ,M/G,6X,9HSPEEU,RPS,/ ,IH ,G 17 .5 , G23 .5,G28.5, 372 

1 G19.5, /,1H0,3X,23HGaS LUnSIANT,J/KG-UEG K, lOX ,8HL, MET ERS , 15X , 373 

2 VHAKEa ,M 2,i4X,5riV,M/S,/,lH , 020. 5 ,G28 .5 , G22. 5 , G20 . 5 ) 374 

51 format IIH ,4X,6HMET£RS,5X,6HKo/i>EC,5X,4HSCMS,6X, 5HCMAX),25X, 375 

1 £ri.MUMtJER,5X, 7HNEwTON0,4X,6MMETeRSJ 376 

52 Format (IHO.lbHMtAN film, meters, 7X,UHPUWER, watts, 8X, 37t 

I 16riT0TAL HEAT, WATTS, 6X,12HDcL( T) ,OEG K,8X,iOHTORQU£, N-M) 378 

5d FORMAT (IHu.llHMEAN FILM = ,01 4. 5 , 2X ,6 HMET ERS, / , IHO , 5X , 8 HX , MET ERS, 379 

1 4X,ariP/PIMlN), 7X,6HP,N/M2,6X, llHUIAVl, M/SEC, 5X,7HMACH NO) 380 

54 FORMAT IlHU,12Aoi 381 

55 FUKMAI I IH , 3X, 6H I NCH £d,6 X ,dHLB/M IN ,5 X ,4HSCFM, 6X, 5H( M AX ) , 25X, 382 

1 EHN JMdcR , 8X , 2riLd , oX , oril NCH ES ) 383 

■Jb format (lM+,99X,3ridAK,7X,3H6ARl 384 

57 FORMAT I IriiJ 385 

C 386 

ENU 387 


IBFIC PRESK 388 

389 

PRESSURE FUNCTION 390 

3<9l 

FUNCTION PRESS! X) 392 

COMMON/ INTGrL/PREF ,RK l , £A , AA, S 393 

EXTERNAL ROT 394 

C 395 

R 1 = KR 1 396 

R = 0 397 

uO = 0. 398 

iF iX.NE.O.J gu=SlMPS2lRl ,u. ,X,RUT ,KT 399 

0 = EXP K(-( X+R I ) #»2 J * (AA* ( UC/SJ +EA) 400 

IF (r.NE.OJ write 16, lU) X 401 

lu FORMAT! IHO, 7HFUR X =, G 1 5. 5 ,3X ,32H I NTEGRAT I ON IN P**2 IS INCORRECT) 402 

Y = S(JRT!U) 403 

PRESS = Y 404 

Return 405 

END 406 
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iibfro expp ^0 7 

c ^08 

C tXPUi'Jt.MTIAL 409 

c 410 

PLfMCIIUN tXPK(X) 411 

uLMMON/tCUN/liCuNil 412 

YY = X 413 

Y = cXP (8CUNir*=YY) 414 

tXPX = Y 415 

JtbUG YY,eCUN8r,Y 416 

KETUKJM 417 

tNU 418 


$IbfTC cXPX 419 

f 420 


C EXT8KNAL FUNCriUN PUK IiMTtGKAL t XP ( E ( X+R 1 1 **2 ) / ( X +R1) OX 421 

C 422 

FEtJOriUiM RUr(Kl,X) 423 

YY = X+Kl 424 

Y = £XPX( YY#*2> /YY 425 

NUI = Y 426 

RETURN 427 

ENU 428 


ilttPTL EPOX 

c 

L external EUNLTIUN fur integral {P-P(REFn UX 

c 

FLnCTUN PX(Y) 

INI LjKL / HKbl“ fAtBfUfU 

yy=/ 

h»x = yy)-HKtF 

K tzT uKN 

hNU 


429 

430 

431 

432 

433 

434 

435 

436 

437 
4 38 


FPXJX 

Q 440 

C EXTERNAL FUNCTION FUR INTEGRAL (P-P(REF))X OX 441 

C 442 

FLNCTIUN PxXtY) 443 

CLNNUN/ INTGRL/PrEF ,A ,D tC,0 444 

yy=/ 445 

PXX = ( PREiil YY) -PREF) *YY 446 

RETURN 447 

cNO 448 
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SMPSKl 449 

450 

NLMEkIC,AL INTcGkAnUN BY SI kPSONS RULE 451 

452 

FLiNiCriUN SI MPSl I XMi N, XMAX,EUf^Cl ,KERJ 453 

JIMcimSIUIm V(2U'J) ,H(200) ,A(20U) ,B(200i ,C { 20U ) , P( 200 ) , E ( 200 ) , NE( 200 ) 454 

EOuIVALENLE {£ ,NE I , ( TE ST.NIESn 455 

T=3.UE- 5 456 

vm=Xi«iIN 457 

HU»=0. 5*(XMAX-XMIN) 458 

A (IJ =fUICl(Xi'1INI 459 

b(iJ=FUNC lUMIN+HUI) 460 

0 ( i)=FUN0 i(XMAX) 461 

P U) =H( U*(A< i)+4.0*oa I + UU » 462 

Em=P(U 463 

AiMS = PU) 464 

'1 = 1 465 

FKAC = 2. C* T 466 

1 FKAt =0. 5 *FkAC 467 

2 rEST=AbSlFi<AL*ANS) 468 

K>=N • 469 

J 00 I 1=1, K, 470 

4 IF (NltSI-lAD j(N£ (11 ) ) 5,5,7 471 

5 N = K+i 472 

V U'll=V( 1 l+n(I 1 473 

HliM 1 = 0. S’f'H ( 1 1 474 

A (N )=Bl II 475 

B {(M 1=FjNC i( V(N)+H IN) 1 476 

C (N) =C( II 477 

P (N)=H(N)*(A(N)+4.0*B(N) + C(N1 ) 478 

0=H(I1 479 

H 11 ) =H( N) 480 

B( i )=FUN0 11 VII ) + rt (I 1 1 481 

Lll)=A(N) 482 

P (1 1 =h( 11*1 A1 I (I H-Ll 11 1 483 

W =P ( I 1+P( Nl-0 484 

ANb=aNS+w 485 

Ed 1=0 486 

E(N)=U 487 

o IF (N-20U1 1,13,13 488 

7 tONriNUE 489 

tJ IF IN-X) 9, 9,2 490 

V 0 = 0.0 491 

10 JC II I =1,N 492 

1 1 J =U<-E (1 1 49 3 

12 IF lAbS(0)-T*AbS(ANS) 1 14,14,1 494 

13 X.cK=RtK + l 495 

14 ANS=0.0 496 

15 00 lb 1=1, N 497 

lo ANb=ANS + Pd 1 498 

SIMP S1=1ANS+0/30. 0) /3. 0 499 

17 HETJKN 500 

ENU 501 
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SHPbRZ 50 2 

503 

iMCHtKlCAL ilM T t^KA Ti Un bY SI MP^iONS RULE 504 

505 

F UMo f lUN SI MP S2 I Y t XMl iM , XMAX ,F UNC2 , kEH ) 506 

u IMcNSlUN V( 200) ,H(2U0J ,A(200J ,b(200i ,CC200), P(200),E( 200),NE(200> 507 

toui VALcNCt (t,Nt j ,( rtST,NTEST) 508 

T = b.uE- 5 509 

V(1)=XMIN 510 

dtU=0, 5*(XMX-XdiNJ 511 

A(iJ=FU.MC2(Y,XMlN) 512 

6 { i i =FUIMC 2 { V , XMI N+H ( U i 513 

C I i) =FUN0 2( Y, XHAX) 514 

P I l)=H( 1)*(A( U<'4-0i^b{l i^C( U ) 515 

E(li=P(l) 516 

ANS=P(1) 517 

iM=i 518 

PKAC=2.0«T 519 

1 FKAO =0. S^'FKaC 520 

2 TE:>T=AbS(FRAC=«^ANS) 521 

522 

3 OU 7 I^i.K 523 

4 If ( wTE ST-IAbS(iMEm ) ) 5,5,7 524 

5 d = N+1 525 

V(iM)=V( U+HU ) 526 

HiN)=0. 527 

A(Ni=b(U 528 

b (N)=HUiNiC2(Y, V(N)+H(i\i ) 529 

C(Ni =C1 I) 530 

P (A(W)+4.0<'b(N)'*-C(N) ) 531 

g=P(I) 532 

d(n=H((M) 533 

b U J =FUNC 2 ( Y , V ( I } +H ( n ) 534 

Cm=A<N) 535 

P m =H( I)^(Am<-4.0*b(l )^-C(U ) 536 

g^P( 1 ) + PUM)-w 537 

ANS=ANS+U 538 

E(I)=U 539 

t ( d ) =g 540 

b IF (.'4-200) 7, lb, 13 541 

7 CUiNillNUE 542 

8 IF ( N-K. ) 9,9,2 543 

90= U.O 544 

10 OU 111 =1 ,N 545 

11 0=g+E(i ) 546 

12 IF ( AbS(g)-T^AbS( ANS) } 14,14,1 547 

13 KEk^KEK'*-/: 548 

14 ANS=0*0 549 

15 DO lo 1=1, N 550 

16 ANb = AUS+P( 1 ) 551 

blMPS2=(ANS+g/3 0. 0) /3.0 552 

17 REIUkN 553 

END 554 
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$IBFTl rkang 

:>LtiKUUf TC AKKANG6 AKkAYS TU b£ 

StdkOUTfNt AKKNG ( X » Y , XR , YP , N , 1 ) 
UiMtilNlGlUN X(^d) ,Yi2ii)J ,XP(^5) , YP(25) 
I = 1 
Ou 100 

If IXIJKLF.U. j GO ro 100 
XPU ) = UJ} 

YP(i ) = X( JJ 
1 = I + i 

100 CUN J iNuE 

1 = i-i 

CALL SORT XY (XP,YP,li 
I A = 1/2 
00 iOi J=i,U 
IJ = 

r = YPijA 
ypu ) = YP( I j) 

YP(iJ) = T 
T = XPl J) 

XPlJj = XPdJi 
XP( 1 J } = 1 

101 CONT INUE 
kETUKN 
tNU 


555 

556 

PLOrTEO 557 

558 

559 

560 

561 

562 

563 

564 

565 

566 

567 

568 

569 

570 

571 

572 

573 

574 

575 

576 

577 

578 

579 

580 

581 
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APPENDIX C 


SAMPLE PROBLEM 

An example of the use of the computer program will now be given with the following 
conditions: Air at 45. 0 psia is to be sealed from ambient air at 15. 0 psia using a seal 
operating in the externally pressurized mode. The mean temperature is 100° F. The 
sealing dam outside diameter is 6. 630 inches, inside diameter is 6. 530 inches, and the 
design speed is 2795 rpm. 

It is desired to find a design mean film thickness which is large enough so that power 
dissipation and viscous heating temperature rise are sufficiently low, yet small enough 
so that the mass leakage is tolerable. From our experience, the best method is to try 
mean film thickness inputs of 0. 1 to 1.0 mil in increments of 0. 1 mil. The output de- 
sired is the mass and volume flow rates, sealing dam force due to the pressure drop, 
center of pressure, power loss, and approximate temperature rise due to shearing. A 
check on the validity of the analysis is made by examining the Mach number, Knudsen 
number, and rotational and pressure flow Reynolds numbers. Thus the program input 
will include 


Speed of rotation, N, rpm 

2796 

Molecular weight of gas, m. lbm/(lb-mole) 

28.9660 

Pressure at inner radius or inlet, psia 

15.0 

Pressure at outer radius or exit, P 2 , psia 

45.0 

Absolute or dynamic viscosity, ju, Ibf-sec/ft^ 

3 .O 6 X lO"”^ 

Specific heat at constant pressure, Cp, Btu/lbm-^R 

0.24 

Mean film thickness (increase in increments 

0. 1 to 1.0 

of 0. 1) , hj^, mil 


Temperature, T, 

100 

Inner radius, R^, in. 

3.265 

Outer radius, R 2 , in. 

3. 315 

Specific heat ratio, y 

1.4 


The data sheet for this sample problem is shown in table HI. The sample output, 
with all possible output options are shown in both English and SI vmits. Note that the 
analysis is invalid for film thicknesses greater than 0. 5 mil (where the Mach number 
has exceeded l/v^= 0. 845). The pressure and velocity distributions, torque, and plots 
of pressure and power versus mean film thickness are included for a more detailed ex- 
amination of this problem. This problem ran in approximately 0. 1 minute on the Lewis 
computer. 
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COMFRESolbLt i^tAUNG UAM wllri RUTAUON AIMO PARALLEL SURFACES 


SAMPLE PKUiiLtiM 
INPUf DATA - 

P^,P:>1A Pl.PSiA r,OEG F VISCOSITY, LB-SEC/FT2 MOLECULAR WEIGHT CP/CV 

4S.00UL iS.OOOu 100. 0.38600E-06 28.9660 1.40000 

R2, INCHES HI, INCHED LYNCHES RH0,L8-’SEC2/Fr4 RHO (RO T > , LB- SEC2/F T4 NO OF GRID POINTS 

3.3liOC 3.2oi>00 0 0 0 U 

N,KPM VfFI/iEc CP,BTU/LB-0£G R SKIP A SKIP R SKIP T 

2796. OC 0 0.24000 F F F 

BEGIN OUTPUT UATA 

RHU(ROr) ,LB-S£C2/Fr4 KHO U ) , LB-SEC2/FT4 AISOJNO SPEED), FT/SEC SPEED, RPM 

0.i)0i2i£-02 0.22470E-02 1160.08 2796.00 

GAS CUNSlANr,Fr-Ld/LB(M)-DbG K L, INCHES AREA, IN2 V, FT/SEC 

5i.ib22 20. 6717 1.03358 80.2750 


MEAN FILM 

M(UUT) 

U 

MACH 

RE(P) 

RE(R) 

KNUOSEN 

F XC 

XC 

F 

INCHES 

Lb/MIN 

SCFM 

(MAX) 



NUMBER 

LB INCHES 

BAR 

BAR 

O.lOuE-02 







ANALYSIS NOT VALID 



t. 900 £-03 







ANALYSIS NOT VALID 



C.aOOE-03 







ANALYSIS NOT VALID 



C.70UE-03 







ANALYSIS NOT VALID 



C.600E-03 







ANALYSIS NOT VALID 



C.500E-U3 

-0.242 

- 3. lo6 

0.670 

377.006 

43.605 

0.005 

18.125 0.318E-01 

0. 635 

0.585 

C.450E-03 

-0.176 

-2.308 

0.543 

274.838 

39.244 

0.006 

18.125 0.318E-01 

0. 635 

0.585 

C.4CUE-03 

-0. 124 

- 1.021 

0.429 

193.027 

34.884 

0.007 

18.125 0.318E-01 

0.635 

0.5 85 

C.350E-03 

-O.B30E-01 

-1.086 

0.328 

129.313 

30.523 

0.008 

18.125 0.318E-01 

0. 635 

0.585 

C.200E-03 

-0.523E-01 

-0.684 

0.241 

81.433 

26.163 

0. 009 

18.125 0.318E-01 

0. 635 

0.585 

C.250E-03 

-0.302E-01 

-0.396 

0.167 

47.126 

21.802 

O.Oll 

18.125 0.318E-01 

0.635 

0.585 

C.iiCOE-03 

-0. 155E-C1 

-0.203 

0.107 

24.128 

17.442 

0.013 

18.125 0.318E-01 

0. 635 

0.585 

C.150E- 03 

-0.653E-02 

-0.853E-01 

0.060 

10.179 

13.061 

0.018 

18.125 0.318E-01 

0.635 

0.585 

C.lOOE-03 

-0. 194E-02 

-0.253E-01 

0.027 

3.016 

8.721 

0.026 

18.125 0.318E-01 

0.635 

0.585 


EAN F1LM,INCHES 

PUWEKjH.P. 

SHEAR HEAT,BTU/M1N 

Dam, DEG F 

TORQUE, FT-LB 

C.1U0E-U2 




ANALYSIS NOT VALID 

C.900E-03 




ANALYSIS NOT VALID 

C.800E-03 




ANALYSIS NOT VALID 

C.700E-03 




ANALYSIS NOT VALID 

C.60UE-03 




ANALYSIS NOT VALID 

C.5C0E- 03 

0. 77908E-03 

0.33048E-01 

0.56909 

0.9195 lE-02 

C.450E-03 

0.B6564E-03 

0.36720E-01 

0.86739 

0.10217E-01 

C.400b-03 

0.9 7385E-03 

0.4131 lE-01 

1.38939 

0.11494E-01 

C.350E-03 

0.11130E-02 

0.47212E-01 

2-37023 

0.13136E-01 

C.300E-03 

0.129B5E-U2 

0.5508 lE-01 

4.39115 

0.15325E-01 

C.250E-03 

0. 15582E-02 

0.66097E-01 

9.10548 

0.18390E-01 

C.200E-03 

0. 19477E-02 

0.82621E-01 

22.2302 

0.22988E-01 

C.150E-03 

0.2596 9E-02 

0.11016 

70.2583 

0.30650E-01 

C.lUOE-03 

0.38954E-02 

0.16524 

355.683 

0.45976E-01 


GO 

CJl 



00 


MtAlM FILM = u.bUoh-Oj 


Xf iNCriL b 

P/P (Mlfyj) 

PrPbl 

U(AV) ,Fl/ShC 

MACH NO 

0 

L • UOOUC) 

lb. UOOO 

777 . lb? 

0 .66992 

O.SOOt- C2 

1.34366 

20. Ib32 

376.380 

0.49857 

0 . iuufc- Ci 

1 . olb46 

24.2310 

48 1.074 

0.41469 

0 • ib Jc- C i 

1 .84730 

^7. 71 

420 .665 

0.36264 

C.ZUuE- C L 

/.Ub^94 

3 0. 7940 

376.559 

0 .32632 

O.^bOb-Ci 

^ .23^^43 

33.b91 7 

347.030 

0.299 14 

C*3out'~ Cl 

2.4113b 

36.1700 

322.294 

0.27762 

0.3b0E-Ci 

2.b715i 

38.D72 7 

302 .216 

0.260 52 

0.40Uc~CL 

2.722U7 

40. 6310 

285.502 

0.246U 

C.4:>Ut-Cl 

2 . 8c4 bl 

42.V677 

271 .305 

0.23387 

U.bLUfc- CL 

3.0U0U0 

4b. 0000 

259 .052 

0.22331 


MfcAN FILM = U.^bOt-Ci INLHfcS 


Xf INCHLS 

P/P(MlNi 

PfPSI 

U(AV> ,FT/6EC 

MACH NO 

0 

1.00000 

15.0000 

629.497 

Q. 54263 

0.5UOb-C2 

1 • b43o8 

2 0. 1552 

466.488 

0.403 84 

O.iOOt-Cl 

1 • 61546 

24.2319 

369.6/0 

0.33590 

C.150fc-Ci 

1.64 730 

27. 7104 

340 .755 

0.29373 

Q.20Ut-Cl 

2.05294 

30. 7940 

306 .633 

0.264 32 

0.2506-U 

2.23945 

33. 591 7 

281.095 

0.24231 

0.300t-U 

2.4113b 

36.1700 

261.058 

0.22504 

0.350t~Cl 

2.57151 

38.5727 

244.796 

0.21102 

0.40ot-’ C 1 

2.72207 

4 0. 6310 

231.257 

u. 19935 

0.45Ufc- Cl 

2.86451 

42. 9677 

219 .757 

0.18943 

0.5OOF-C1 

3.OCOO0 

45. 0000 

209.832 

0.18088 


MfcAN FiLM = 0.^00£-03 INCHtS 


X, INCHES 

P/P(M1NJ 

PfPSl 

U(AV) , FT/SEC 

MACH NO 

0 

l.OGOOG 

1 5. UO Ou 

497.381 

0.42875 

0.5U0t-C2 

1.34360 

20.1552 

370 .163 

0.31909 

0 . loot- C i 

1.61546 

24. 23 1 9 

307 .688 

0.26540 

0.150E-C1 

1 . 84736 

27. 7104 

269.239 

0.23209 

0.2U0E-C1 

2.05294 

30. 7940 

242.276 

0.20885 

0.250E-C1 

2.23945 

33.5917 

222.099 

0.19145 

0.300fc-Ci 

2.41133 

36.1700 

206.268 

0.17781 

0.350E-C1 

2.57151 

38.5727 

193 .419 

0 . 16673 

0.400t-Cl 

2.72207 

40.6310 

182.722 

0.15751 

0.450t- Cl 

2.86451 

42.9677 

173.635 

0.14968 

0.500t-Cl 

i.OOOOO 

45.0000 

165.794 

0.14292 



MtAN FILM 


U . 3 5\Jt- C3 I iMLrlfc 5 


Mb AM 


MtAM 


Af iNLHt S 

P/P (MliM) 

P,PSI 

OIAVi , FT/SEC 

MACH NO 

0 

i .00000 

15.0000 

380.807 

0.32826 

0.300b-C^ 

i .34350 

20.1552 

283.406 

0.24430 

U- iU0£- L 1 

i . oi545 

24.2319 

235.726 

U- 203 20 

O.llsUb-Li 

i.d473o 

27.7104 

206 .136 

0.17769 

U. 20ut- C i 

2.05^94 

3 0. 7940 

185.494 

0.15990 

0 • 2!> J£— C i 

2.23945 

33. 591 7 

170.045 

0.146 58 

0.3UUL-C1 

2.4ii33 

36.1700 

157.924 

0.13613 

0.330t- U 

2.37i5i 

30-5727 

148.087 

0.12765 

0-4U0l“ C 1 

2. 7^207 

40. b310 

139.896 

0.12059 

U .4b Ob'- C 1 

2.06451 

4 2 . 90 7 7 

132.940 

0.11460 

U . 3UUL- C I 

3.00000 

45. 0000 

126.936 

0.10942 

F iL M = b . 

300E-U3 INCHES 




X, ii\CH£S 

P/P(MIN) 

P,PSi 

U(AVi ,FT/SbC 

MACH NO 

0 

1 .00000 

15.0000 

279.777 

0.24117 

0.5u0t- G2 

i .3430d 

20. 1552 

208.217 

0.17949 

U,lUOb-Ci 

i.oi34o 

24,2319 

173.187 

0.14929 

C.lbO£-Ci 

i .04736 

27.7104 

151.447 

0.13055 

G.2LiUt-Ci 

2.05294 

30. 7940 

136.281 

0.117 48 

0.2bUb“Ll 

2.23945 

33.5917 

124.931 

0.10769 

0,300E“CI 

2.41i33 

36.1700 

116.026 

0.10002 

0.350b- Ci 

2.5715i 

38. 5727 

108.798 

0.93785E-01 

0.40ot- Ci 

2.72207 

4 0. 8310 

102. 781 

0.88598E-01 

C.450b-Ci 

2.0645i 

42. 9677 

97.6699 

0.84193E-01 

C.bOOb- Ci 

3.00000 

45.0000 

93.2588 

0. 80390 E-Oi 

FiLM = 0. 

250E-03 INCHES 




X, iiMOHbS 

P/P IMINi 

P,PS1 

U(AVi , FT/SEC 

MACH NO 

0 

i.oooou 

Id. 0000 

194.289 

0.16748 

0.5OOE-C2 

i. 34360 

20. 1552 

144.5 95 

0.12464 

0, iOUb- Ci 

i.6i546 

24.2319 

120.269 

0.10367 

C.i50b-Ci 

1.04736 

2 7. 7104 

105.171 

0. 90659 E“Ol 

0.200£-Ci 

2.05294 

30.7940 

94.6397 

0.81581E-01 

C.250E-Ci 

2.23945 

33.5917 

86.7576 

0.747 86 E-01 

0.300b-Ci 

2.4U33 

36.1700 

80.5735 

0.69455E-01 

0.35Ub- Ci 

2.57151 

38.5727 

75.5544 

0. 65129 E-01 

0.4UUb- Li 

2.72207 

40.8310 

71.3756 

0.61527E-01 

C.450t-Ci 

2« 86451 

42.9677 

67.8263 

0.58467E-01 

0.500£-Ci 

3.00000 

45. 0000 

64.7631 

0.55827E-01 


00 

-Cl 



00 

00 



UiAV) , FT/SEC 
124.345 
92.5409 
76.97 19 
67.3U96 
&0.569.4 
55.5249 
51.5670 
H8.354B 
45.6804 
43.4088 
41.4484 


U(AV) ,FT/StC 
69.9441 
52.0542 
43.2967 
37.8617 
34.0703 
31.2327 
29.0064 
27.1996 
25.6952 
24.4175 
23.3147 


U(AV) , FT/SEC 
31.0863 
23.1352 
19.2430 
10.8274 
15.1424 
13.8812 
12.8918 
12.0887 
U.4201 
10.8522 
10.3621 


MACH NO 

0.10719 

0.79771E-01 

0.66351E~01 

0.58022E~01 

0.52212E-01 

0.47863E-01 

0.44451E-01 

0.41682E-01 

0.393776-01 

0.37419E-01 

0.35729E-01 


MACH NO 
0.60293E-01 
0.44871E-01 
0.37322E-01 
0.32637E-01 
0.29369 E-01 
0.269 23 £-01 
0.25004 E-01 
0.23446 £-01 
0.22150 £-01 
0.21048E-01 
0.20098E-01 


MACH NO 

0.26797E-01 

0.19943E-01 

0.16588E-01 

0.14505E-01 

0.13053E-01 

0.L1966E-01 

0.1U13E-01 

0.10421E-01 

0.98443E-02 

0.9 3547E-02 

0.89322E-02 



PLOT OF POriER VS H(MEAN) 


U. OOUiOU 

U. OU^OOl 

* 

i 

i 

I 

I 

I 

I 

1 

I 

0. C030U1 
I 
1 
1 
1 
I 
I 
1 
L 
I 

0. UQ2001 
i 
1 
1 
I 
I 
I 
1 
1 
I 

O.OOIOOl 

I 

1 

1 

i 

I 

1 

i 

I 

I 

a. i — 

O.UOQIOO 


0. 000150 0 .000 200 0.0 00 250 0.000300 0.000350 0.000400 0.000450 0.000500 0.000550 


0. 000150 0.000200 0.000250 0.000300 0.000350 0.000400 0.000450 0.000500 0.000550 

POMER ON VERTICAL SCALE - H(H£AN) ON HORIZONTAL SCALE 

POHER IN HORSE POWER - TO CONVERT TO WATTS, MULTIPLY BY 745.7 
H(MEAN) IN INCHES - TO CONVERT TO METERS, MULTIPLY BY 2.54E-2 


0.000600 


0.000600 


FOR SHEAR HEAT IN BTU/MIN, MULTIPLY POWER BY 42 .42 



O 


PLOT UF P/P (MlN ) Vb X FOR H( ME AN) = 0.5u000t-J3 

0. U.OOdU O.OlOO 0.0iS)0 

I 

I 

I 

L 

1 

1 

I 

1 

I 

Z.^Oi 

i 

I 

i 

i 

I 

I 

1 

1 

1 

2,00i 

1 

1 

1 ♦ 

1 

i 

I 

1 

I * 

1 

1.501 

1 

1 

1 * 

1 

1 

1 

1 

1 

1 

1 . 00 * 

1 

1 

1 

1 

1 

1 

1 

1 

1 



0 . 0 . 005 0 0.0 10 0 0 . 0 1 50 

P/PCMINJ ON VtKIlCAL SCALE - X ON HORIZONTAL SCALE 


0.0200 0.0250 0.0300 


* 


* 


* 


0*0200 0.0250 0.0300 


0.0350 0.0400 0.0450 


* 


♦ 


* 


0*0350 0.0400 0.0450 


0.0500 


0.0500 


P/PIMIN) IS DIMENSIONLESS 

X IN INCHES - TO CONVERT TO METERS, MULTIPLY BY 2.54E'-2 



C0MPkE:>SI6LE StALXNG JAM WITH KU TATI UN ANU PARALLEL SURFACES 


sample PRUdLEM 


P2.N/M2 

Pi,N/M2 

T* JEG 

K 

VISCOSITY 

,N-S/M2 

MOLECULAR WEIGHT 


CP/CV 

0 

C6 0.10J42E 06 

311 

• 

0. 18482E-04 

28.9660 


1.40000 

K2,MtTER 

S Ri, METERS 

L , ME I E RS 

KHQ,KG/M3 

RH0(R0T),KG/H3 

NO 

OF GRID POINTS 

0.8^2UiE- 

Cl 0.829ilt-0l 

0 


0 


0 



11 

N , RP S 

V,M/S 

LPf J/KG 

.-DEG K 


SKIP A 

SKIP R 

SKIP T 



A6*6000 

0 

1004. 

78 


F 

F 

F 



BEGIN UUIPUT 

UATA 









KHU(KUT),RG/M3 RHUU) , 

RG/M3 

AISGOND SPEED), M/S 

SPEED, KPS 




Z.GU26i 1.16217 



333.591 


46.6000 




GAS LUN6TAN r f J/KG— JEG K 

L r ME T£K6 

AREA, 

M2 

V.M/S 




287 

• 086 

0.52506 


0.66683E-03 

24.4678 



MEAN FILM 

M(JUT) g 

malh 

RE(P| 

RECRI 

KNUOSEN 

F 

XC 

XC 

F 

METERS 

KG/SEC SCMS 

IMAX) 



NUMBER 

NEWTONS 

METERS 

BAR 

BAR 

0.25<»E-04 






ANALYSIS 

NOT VALID 



C.229E-04 






ANALYSIS 

NOT VALID 



0.203E-04 






ANALYSIS 

NOT VALID 



0. I7ttE“U4 






ANALYSIS 

NOT VALID 



0.152E-U4 






ANALYSIS 

NOT VALID 



0.127E-04 

-0.183E-02 -0. 149E-02 

0.6 70 

377.006 

43.605 

0.005 

80.623 

0.807E-03 

0.635 

0.585 

0.114E-04 

-U.I33E-C2 -U.1U9E-02 

0.343 

274.838 

39.244 

0.006 

80.623 

0.807E-C3 

0.635 

0.585 

0.102E-U4 

-0.937E-03 -0. 7oSE-03 

0.4Z9 

193.027 

34.884 

0.007 

80.623 

0.807E-03 

0. 635 

0.585 

C.fa89E-0S 

-U.627E-03 -U.512E-03 

0.328 

129.313 

30.523 

0.008 

80.623 

0.8C7E-03 

0.635 

0.585 

0.762E-U5 

-0.393E-C3 -U.323E-03 

0.241 

81.433 

26.163 

0.009 

80.623 

0.807E-03 

0. 635 

0.585 

0.63SE-06 

-0.229E-U3 -0.187E-03 

0.167 

47.126 

21.802 

O.Oll 

80.623 

0.807E-03 

0. 635 

0.585 

0.5C8E-05 

-O.lirE-03 -0.936E-04 

0.107 

24. 128 

17.442 

0.013 

80.623 

0.8C7E-03 

0. 635 

0.585 

0.28 IE- US 

-0.49*»E-C4 -0.40:>E-U4 

0.060 

10.179 

13.081 

0.0L8 

80.623 

0.807E- 03 

0.635 

0.585 

C.254E-0S 

-0.1ho£-C 4 -0.120t-04 

0.027 

3.016 

8.721 

0.026 

80.623 

0.807E-03 

0.635 

0.585 

MEAN FlLMfMclEKS POwEK , wATTS 

ruTAL HEAT 

.WATTS 

DEL(T),DEG 

K TORQUE, N-M 



0.254E-04 






ANALYSIS 

NOT VALID 



0.229E-04 






ANALYSIS 

NOT VALID 



C.2C3E-04 






ANALYSIS 

NOT VALID 



C. 178E-04 






ANALYSIS 

NOT VALID 



0.152E-U4 






ANALYSIS 

NOT VALID 



0.X27E-U4 

0. 58096 


0.38158 


0.31616 

0 

.12467E-01 



0.114E-04 

0.64551 


0.64620 


0.48188 

0 

.13852E-01 



0.102E-U4 

0. 72 62 0 


0.72698 


0.77188 

0 

.15584E-01 



C.889E-US 

0.82 994 


0.83083 


1.31680 

0 

.17810E-01 



(i.762t-u5 

0 . 96 826 


0.96930 


2.43953 

0 

.20778E-01 



0.63SE-U5 

1. 16191 


1.16310 


3.05 860 

0 

.24934E-01 



0.SG8E-U5 

1.45239 


1.43395 


12.3501 

0 

.31167E-01 



0.381E-US 

1.93632 


1.93860 


39.0324 

0 

.41556E-01 



0.254E-05 

2.90479 


2.90790 


197.602 

0 

.62334E-01 





4 ^ 

to 


FILM = o.i^/CCfc-o4 Mi£TEKb 


XfMtThK S 

P /P lMl^4) 

P,N/M2 


UtAV) fM/66C 

MACH NO 

0 

I • OUuOO 

0. i0342E 

Oo 

236.077 

0.66992 

Q.LcJt- 

i « 3^^6d 

0. 1 309 7t 

06 

176.290 

0.49857 

0.2t>4b- Cj 

1.015^6 

0.16707E 

06 

14o .631 

0.41469 

0.3b iL- L3 

i . d^ iso 

0. 191 06 E 

06 

128.225 

0.36264 

0.5C3b- CJ 

2.0b294 

0.21232E 

06 

U5.305 

0.32632 

0.63:>t- C3 

2.2394b 

0.23161E 

06 

105.775 

0 .299 14 

0.76^t- C3 

2.^U33 

0.24930E 

06 

90.2352 

0.27782 

O.dbyt- C3 

2.:>7ibl 

0.26595E 

06 

92.1159 

0.26052 

U.i02c- C2 

2. 722J7 

0.2dl52E 

06 

07.0211 

0.24611 

O.l i4t-C^ 

2.o64!)l 

0.29625fc 

06 

02 .6930 

0.23387 

G.i^7fc-L2 

3.OUU0C 

0.310266 

06 

70.9592 

0.22331 

FILN = U. 

ii^dUE-04 

ME TEk S 




X,WhI tr< S 

P /P (MIN) 

P,N/M2 


U(AV) f M/S6C 

MACH NO 

0 

I .uoouo 

U.10342E 

06 

191.071 

0.54263 

0.127E-C3 

i. 34308 

0.130976 

06 

142.795 

0.40384 

0.23<*t-C3 

1 • blb4fa 

0.Io7076 

Ob 

110.771 

0.33590 

0.3dib-C3 

i. 04730 

0.1910£>£ 

06 

103.862 

0.293 73 

0.50de- C3 

2.05294 

0.212326 

06 

93.4616 

0.264 32 

0.63i)C- C3 

2.23945 

0.231616 

06 

05.6776 

0.24231 

0. /6^fc- L3 

2.4U33 

0. 249306 

06 

79.5705 

0.22504 

0.8d9t- C3 

2.57ibl 

0.265956 

06 

74.6139 

0.21102 

0. i02t- C2 

2. 72207 

0.20152 6 

06 

70.9071 

0.19935 

0.ll4t-C2 

2.36451 

0.296256 

Oo 

66.9820 

0.18943 

0.i2 7t- C2 

3.00UOO 

0.3102o£ 

06 

63.9569 

0.18088 

FiLh = J. 

lJi60E-U4 

MEIERS 




XrMhTFK S 

P/P (MIN) 

P,N/M2 


U(AV), M/SEC 

MACH NO 

0 

1 .00000 

0. 10342 6 

Oo 

151.602 

0.42875 

0.127t-C3 

1.34308 

0.138976 

06 

112.026 

0.31909 

0.2b4E-C3 

1 . 61540 

0.1o707E 

06 

93.8441 

0.26540 

0.381E-C3 

1.04736 

0.19106 6 

06 

82.0639 

0.23209 

0.50dE-C3 

2.05294 

0.21232E 

06 

73.0462 

0.20885 

0.<>3^E“C3 

2.23945 

0.2316LE 

06 

67.6959 

0.19145 

0.762E-C3 

2.41133 

0.249306 

06 

62.8705 

0.17781 

0.d89E~C3 

2.57151 

0.26595E 

06 

58.9542 

0.16673 

0. I02t- C2 

2.72207 

0.281 52 E 

06 

55.69 35 

0.15751 

0.114E- C2 

2.06451 

0.29625E 

06 

52.9241 

0.14968 

0.127E-C2 

3.00000 

0.31026E 

06 

50.5339 

0*14292 



MtAN FiL.-l 


U- bdSCCt-U5 


ME TEKb 


X,ME TEK S 

i* /P (Ml M) 

P,N/M2 


U(AVi, M/SEC 

MACH NO 

0 

i .00000 

0. 10342 E 

06 

116.070 

0 .32826 

G.I2 7E-C3 

i . i4i6 d 

0.li897E 

Oo 

Bo. 3823 

0.24430 

0.21?^E“ C3 

i . o i i 4t> 

0. i 6 7 07E 

06 

71.6494 

0.20 3 20 

0- do it- Ci 

I. d47io 

0.i9iU6E 

06 

62.8302 

0.17769 

u. 5 OaC- C3 

2.US294 

0.21232E 

06 

36.5385 

0.15990 

0 •<? 3l?E- C3 

2 .2i94i> 

0.2il6lE 

06 

31.8297 

0.14658 

0- /o2h- Ci 

2.4iiii 

U.249idt 

06 

46 .1332 

0.13613 

0.B89E- Ci 

2.i)711?i 

0. 2o5 *y3t 

06 

45.1368 

0.12765 

0.1u2t- C2 

2. /2207 

0.261 32E 

06 

42.6404 

0.12059 

O.i Uc-C2 

2 .do4il 

0.29623E 

06 

4o . 520u 

0.11460 

0.127E-L2 

i. 00000 

0. 31026E 

06 

38.6900 

0. 10942 

F iL M = u . 

7O2O0t-Ui 

Mt TtRS 




XiMETEK S 

P/P (MINI 

P,N/M2 


U( AVI fM/S£C 

MACH NO 

0 

i .00000 

0. 10342E 

06 

85.2739 

0.24117 

O.L2?E-Ci 

i .i4ido 

0.i3697E 

06 

63.4645 

0.17949 

0.234E-C3 

1.61346 

0.16707E 

06 

32.7873 

0.14929 

U. JbiE-C3 

L.dH73o 

0.1910oE 

06 

46.1609 

0.13055 

0.5UdE- Ci 

2.03294 

0.21232E 

06 

41.3385 

0.11748 

O.bibE- C:> 

2.2i943 

0.2 3161E 

06 

38.0789 

0.10769 

0.762E-Ci 

2.4liii 

0.24936E 

06 

35.3647 

0.10002 

0.db9t- Ci 

2.37131 

0.26393E 

06 

33.1617 

0.93785E-01 

0,iU2E-C2 

2. 72207 

0.26132E 

06 

31. 3276 

0.885986-01 

0.li4t-C2 

2.86431 

0.29623E 

06 

29.7698 

0.841936-01 

0.127E-C2 

i.oOOUO 

0.31026E 

06 

28.4253 

0. 80390 E-01 


MEAN FILM = U . bOOE- 01? METERS 


X,MtTERS 

P/P (MINJ 

P,N/M2 


U(AV) , M/SEC 

MACH NO 

0 

1.00000 

0.10342L 

06 

39.2194 

0.16748 

0.127E-C3 

1.34368 

0.13897E 

Oo 

44.0726 

0.12464 

0.254E-C3 

1.61346 

0.16707E 

06 

36.6579 

0.10367 

0.381t-C3 

1.84730 

0.19106E 

06 

32.0562 

0. 90659 E-Ol 

0.5C8E-C3 

2.05294 

0.21232E 

06 

28.8462 

0.81581E-01 

0.635t-Ci 

2.23943 

0.231ol£ 

06 

26.4437 

0.74786E-01 

0.762E-C3 

41133 

0.24938E 

06 

24.5588 

0.69455E-0I 

0.889E-C3 

2.37151 

0.26593E 

06 

23.0290 

0.65129 E-01 

0.102E-C2 

2.72207 

0.28152E 

06 

21.7553 

0.61527E-0I 

U.114E-C2 

2.86451 

0.2 9625E 

Ot> 

20.6735 

0.58467E-01 

0.127E-C2 

3.00000 

0.31026E 

06 

19.7398 

0.55827E-01 


CO 
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MhAiN FILM 


U. bOcOLc-Lb MtIbKb 


X#MLTtR b 

r* /P IMI R) 

P, N/M2 


JIAV) , M/SEC 

MACH NU 

u 

1 . UOJUU 

0. 10342 E 

06 

37.9J04 

0.10719 

0 - tZ7t-C:> 

i • 3 ^3od 

0. 1 3o97E 

06 

28 . 20 Ob 

0 .79771E-01 

0 . 2 r>^t- CJ 

I -b Ib-HO 

O.16707L 

06 

2J.46 10 

0.66351E-01 

0 - :>fc Lt- C3 

1 . o4 V JO 

0.191 06 E 

00 

20.3160 

O.58022E-01 

U.^Lbfc- L J 

2,0b.iSi4 

0.21232E 

06 

18 .4o 16 

0.52212E-01 

U .o 3bL- C 3 

2 .23S48 

0.23161E 

Oo 

16.9240 

0.47863E-0L 

u, / o 2 t- L j 

2 .41 1 J3 

0.2493 8 E 

Oo 

13.7176 

0.44451E-01 

u,d89h-L3 

2. 37131 

0.265 93E 

06 

14. 7383 

0 .41682E-01 

0 . L 02 t- 

722U7 

0.281 32E 

06 

13.9234 

0.39377E-01 

O.I 

2.86431 

0.29O23E 

06 

13.2310 

0.37419E-01 

U.U7h-C2 

3 . oOOuO 

0. 31026b 

06 

12.6333 

0-35729E“01 

FILM = 0 . 

3810LE-03 

Ml I tHb 




XtMEFLKb 

P/P (MINI 

P,N/M2 


UIAV), M/SEC 

MACH NO 

0 

1 .UOOOO 

0.10342E 

06 

21.3190 

U.60293E-01 

O.L27tf~L3 

1 .34308 

0.13d97E 

06 

13.8661 

0 .44871E-01 

U .2ij4£- C3 

i . 0 I 34o 

0.16707E 

00 

13. 1968 

0.37322E-01 

U.381L-C3 

i . 84 730 

0.19106E 

06 

11.3402 

0.32637^01 

0.bU8t-C3 

2.05294 

0.21232b 

06 

10.3846 

0.29369E-01 

0.63i>t- C3 

2 .23945 

0.23161E 

06 

9.31974 

0.26923E-01 

0.762E-C3 

2.4U33 

0.249386 

06 

8.84117 

0.25004E-01 

0.dbVt-C3 

2.37131 

0.26595E 

06 

8.29043 

0. 23446 E-Oi 

0 . 102 E-C 2 

2.72207 

0.281 32E 

06 

7.83190 

0. 22150 E-Ol 

U.114E-C2 

2.86451 

U. 29625b 

06 

7.44245 

0.21048E-01 

0.l27t- C2 

3.00GU0 

0 . 31026E 

06 

7.10632 

0.20098E-01 

FILM = 0 . 

234U0E-03 

METERS 




X, METERS 

P/P IMINI 

P,N/M2 


UlAVI , M/SEC 

MACH NO 

0 

1 .00000 

U.10342E 

06 

9.47510 

0.26797E-01 

0.i2?E-C3 

1 .34368 

0.13897E 

06 

7.05161 

0.199436-01 

0,2b4E-C3 

1 • 61540 

0.16707E 

06 

3.86526 

0.16588E-0L 

0.3BIE-C3 

1.84736 

0.19106E 

06 

5.12899 

0.145056-01 

0,50dE-C3 

2.05294 

0.21232E 

06 

4.61539 

0.130536-01 

0.633E-C3 

2-23945 

0.23161E 

06 

4.23099 

0.119 66E-01 

0.762E-C3 

2.41133 

0.24938b 

06 

3.92941 

0.111136-01 

Q.889E-C3 

2.57131 

0.26393E 

06 

3.68464 

0.L0421E-01 

0.i02E-C2 

2.72207 

0.28132E 

06 

3.48083 

0.98443E-02 

O.I I4h- L2 

2.86451 

0.2 962 5E 

06 

3.30775 

0. 935476-02 

0.127E-C2 

3.00000 

0. 31026E 

06 

3.15837 

0.89322E-02 
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